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Introduction 

This  project  is  to  decrease  the  number  of  breast  cancers  that  are  missed  in  conventional  mammography  by 
combining  two  methods  developed  to  increase  the  sensitivity  of  breast  cancer  imaging:  digital  breast 
tomosynthess  (DBT)  and  temporal  subtraction.  A  major  limitation  of  traditional  mammography  is  the 
confluence  of  overlying  structures  resulting  from  a  two-dimensional  (2D)  projection  of  three-dimensional  (3D) 
anatomy.  DBT  is  an  exciting  new  modality  for  breast  imaging  which  strives  to  overcome  the  obstacles 
presented  in  conventional  2D  mammography  by  taking  multiple  projections  over  a  fixed  angle  and 
reconstructing  volumetric  data  which  serves  to  isolate  overlying  anatomy  from  in-plane  structures.  Such  section 
images  serve  to  amplify  the  conspicuity  of  lesions,  particularly  in  dense  breasts.  DBT  has  been  shown  in  several 
studies  to  increase  the  sensitivity  of  mammography  for  breast  cancer  detection.  One  of  the  primary  methods 
used  by  radiologists  to  detect  developing  tumors  is  comparative  analysis.  This  is  the  process  of  visually 
analyzing  the  temporal  change  between  a  current  and  prior  mammogram,  utilizing  perception  to  perform  the 
requisite  geometric  transformation  between  images.  Temporal  subtraction  is  not  dissimilar  from  comparative 
analysis.  It  automates  the  process  by  using  two  images  taken  sequentially  and  subtracting  them  in  order  to  find 
temporal  discrepancies.  The  purpose  of  this  project  is  to  detennine  the  feasibility  of  using  temporal  subtraction 
on  DBT  phantom  images  to  allow  for  easier  and  earlier  detection  of  breast  cancer  than  with  either  technique 
alone.  A  computer  simulated  phantom  will  be  developed  to  generate  tomosynthesis  data  sets  which  can  be  used 
for  the  development  and  evaluation  of  an  automated  registration  technique  which  will  then  be  applied  to 
physical  phantom  tomosynthesis  images  in  a  graphical  user  interface  for  the  purpose  of  performing  an  observer 
study  to  asses  the  developed  technique. 


Body 

Task  1.  To  generate  tomosynthesis  datasets  of  simulated  and  physical  breast  phantom,  Months  1-12: 

la.  Develop  a  realistic  computer  simulated  breast  phantom  simulation  and  generate  up  to  50  simulated 
tomosynthesis  projection  data  with  the  phantom  undergoing  simulated  tissue  deformation. 

After  investigating  different  types  of  computer  simulated  breast  phantoms,  the  investigator  believes  that  a 
combination  of  a  mathematical  and  voxelized  breast  phantom  would  be  the  best  path  to  pursue  for  the  computer 
simulated  phantom  development.  The  investigator’s  colleagues  have  access  to  high  resolution  breast  CT 
datasets  that  may  be  used  in  the  phantom  development.  After  seeking  the  proper  institutional  approval  (e.g. 
IRB)  the  investigator  plans  to  segment  the  breast  CT  data  with  an  automated  segmentation  algorithm.  The  lab 
of  the  investigator  has  a  technique  which  enables  the  combination  of  multiple  CT  datasets  into  a  computer 
phantom  that  is  flexible  in  its  size  and  shape  called  the  four  dimensional  (4D)  non-unifonn  rational  b-splines 
(NURBS)  based  Cardiac-Torso  (NCAT).  The  segmented  data  will  be  utilized  to  create  a  detailed  3D  computer 
generated  breast  phantom  based  on  empirical  data  using  a  combination  of  non-unifonn  rational  b-splines 
(NURBS)  and  subdivision  surfaces  (SD).  The  phantom  will  be  applicable  to  many  different  types  of  breast 
imaging  research  and  will  be  adjustable  in  size  and  deformable  using  finite-element  methods.  For  the  purpose  of 
this  research  project,  the  3D  computer  simulated  phantom  will  provide  a  realistic  and  flexible  breast  phantom 
for  the  purpose  of  registration  method  development  by  incorporating  simulated  compression  and  the  generation 
of  simulated  tomosynthesis  breast  acquisitions. 

Issues:  The  investigator  needs  to  acquire  institutional  approval  to  work  with  human  subject  data  before 
continuing  with  this  part  of  the  project.  After  approval,  the  investigator  plans  to  use  some  form  of  adaptive 
region  growing  in  addition  to  filtering  and  thresholding  in  order  to  overcome  any  noise  issues  and  make  the 
algorithm  more  robust. 
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lb.  Acquire  up  to  30  tomosynthesis  projections  of  a  compressible  and  deformable  physical  phantom  with 
physically  simulated  anatomy  and  under  different  simulated  temporal  discrepancies. 


The  investigator  used  a  compressible  and  deformable  physical  breast  phantom  shown  in  Figure  1.  Inside  of  the 
breast  phantom  the  investigator  inserted  a  number  of  objects  to  simulate  breast  parenchyma:  acrylic  sponge  and 
yam  (shown  in  Figure  2)  to  simulate  fibroglandular  tissue,  differently  sized  beans  and  acrylic  spheres  to 
simulate  nodules  (shown  in  Figure  3),  small  pieces  of  egg  shell  to  simulate  calcifications,  and  mineral  oil  to 
simulate  adipose  tissue.  The  investigator  took  multiple  tomosynthesis  exams  using  32  kVp,  125  mAs,  and  with 
the  breast  compressed  to  ~7cm  with  a  slight  rotation  in-between  image  acquisitions.  There  were  three  scans 
taken  in  total:  one  in  the  original  orientation  and  original  number  of  simulated  lesions;  the  second  with  the 
breast  slightly  rotated  with  the  original  number  of  simulated  lesions;  and  the  third  with  the  breast  rotated  in  a 
different  way  and  without  some  of  the  simulated  lesions.  Example  projections  from  the  compressed  phantom 
are  shown  in  Figure  4. 


Issues:  The  materials  used  for  the  physical  phantom  were  not  ideal  and  the  tomosynthesis  images  reconstructed 
from  the  acquisition  do  not  look  realistic  enough  to  use  for  the  observer  study.  Further  research  must  be  done  to 
match  up  materials  with  acquisition  parameters  to  create  more  realistic  images. 
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Figure  3:  Projection  of  Physical 
Phantom 


The  lab  of  the  investigator  has  access  to  a  set  of  DBT  human  subject  images  which  contain  pairs  of  temporally 
sequential  data.  In  addition  to  pursuing  the  physical  phantom  pathway,  the  investigator  plans  to  acquire 
institutional  approval  (e.g.  IRB)  to  work  with  the  human  subject  data.  This  will  give  the  investigator  human 
data  to  work  with  for  this  project  and  perhaps  lend  the  technique  developed  a  more  compelling  clinical  and 
practical  implementation. 

Issues:  IRB  approval  must  be  acquired  prior  to  using  the  temporally  sequential  human  subject  DBT  images. 

The  investigator  has  investigated  a  simulated  dose  reduction  technique  for  use  on  tomosynthesis  images  using 
anthropomorphic  chest  phantom  images.  This  research  was  performed  on  chest  data  in  order  to  continue  with 
prior  research  and  validate  the  technique  on  images  the  investigator  had  already  acquired.  We  experimentally 
determined  the  NPS  of  the  tomosynthesis  acquisition  system  and  utilize  it  to  filter  an  image  of  random  noise. 
After  some  further  modifications  to  adjust  the  noise  variance,  this  resultant  noise  image  is  added  to  the  original 
image  and  the  procedure  culminates  in  an  image  which  simulates  an  image  acquired  at  a  reduced  exposure.  This 
technique  is  easily  applicable  to  the  breast  tomosynthesis  data  that  the  investigator  will  acquire  in  the  course  of 
this  project  and  can  be  used  to  additionally  evaluate  the  dose  of  breast  tomosynthesis  images  on  lesion 
detectability.  Please  see  Appendix  1  for  further  details. 

Issues:  The  dose  reduction  technique  that  was  developed  for  the  chest  images  must  be  altered  to  be  used  on 
breast  tomosynthesis  images. 

lc.  Utilize  up  to  3  different  tomosynthesis  reconstruction  algorithms  (Filtered  Back  Projection,  Matrix  Inversion 
Tomosynthesis,  and  Gaussian  Frequency  Blending)  to  create  tomosynthesis  data  sets  of  the  simulated  and 
physical  phantom. 

The  DBT  images  were  reconstructed  using  the  Filtered  Back  Projection  algorithm.  As  mentioned  previously 
the  images  do  not  appear  to  realistically  resemble  actual  breast  tomosynthesis  images  from  human  subjects. 
Please  see  Figure  5  for  examples  of  tomosynthesis  slice  through  the  reconstructed  volume. 

Issues:  None. 
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Figure  4:  Tomosynthesis  slice 
of  physical  phantom 


Key  Research  Accomplishments 

•  Identification  of  simulation  model  to  use  combining  real  human  data  and  geometric  primitives 

•  Physical  breast  phantom  tomosynthesis  images  have  been  acquired  and  processed 

o  New  materials  must  be  identified  for  use  in  physical  breast  phantom  for  more  realistic  images 

•  Tomosynthesis  simulated  dose  reduction  technique  has  been  developed 

o  Needs  to  be  applied  to  breast  images 


Reportable  Outcomes 

1.  C.  M.  Li  and  J.  T.  Dobbins,  III,  “Methodology  for  Detennining  Dose  Reduction  for  Chest  Tomosynthesis,” 
SPIE  Medical  Imaging  2007:  Physics  of  Medical  Imaging  6510,  (2007). 

Conclusions 

The  identification  of  the  proper  computer  simulation  method  is  an  important  milestone  and  will  provide  a 
research  tool,  not  only  for  the  scope  of  this  project,  but  also  for  use  by  other  institutions  who  are  working  on 
breast  imaging.  Although  further  work  must  be  performed  towards  choosing  the  correct  materials  for  the 
physical  breast  phantom,  it  was  important  to  show  that  images  could  be  acquired  and  they  can  still  be  used  for 
the  development  of  the  registration  algorithm  while  the  computer  simulated  phantom  is  being  created.  DBT  has 
been  shown  to  increase  the  specificity  and  sensitivity  of  lesion  detection;  however  dose  reduction  is  an 
important  addition  for  clinical  implementation  of  DBT  to  constrain  the  radiative  dose  to  the  patient  during 
screening.  The  methodology  developed  during  this  project  for  dose  reduction  will  provide  a  good  simulation 
technique  for  evaluating  the  dose  limitations  for  DBT. 
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Digital  tomosynthesis  is  an  imaging  technique  that  reconstructs  tomographic  planes  in  an  object  from  a  set  of  projection 
images  taken  over  a  fixed  angle1.  Preliminary  results  show  that  this  technique  increases  the  detectability  of  lung 
nodules2.  Current  settings  acquire  images  with  approximately  the  same  exposure  as  a  screen-film  lateral.  However, 
due  to  the  increased  detectability  of  lung  nodules  from  the  removal  of  overlying  structures,  patient  dose  may  be  reduced 
while  still  maintaining  increased  sensitivity  and  specificity  over  conventional  chest  radiographs.  This  study  describes  a 
simulation  method  that  provides  realistic  reduced  dose  images  by  adding  noise  to  digital  chest  tomosynthesis  images  in 
order  to  simulate  lower  exposure  settings  for  the  purpose  of  dose  optimization.  Tomosynthesis  projections  of  human 
subjects  were  taken  at  dose  levels  which  were  specified  based  on  either  patient  thickness  or  a  photo-timed  digital  chest 
radiograph  acquired  prior  to  tomosynthesis  acquisition.  For  the  purposes  of  this  study,  subtle  nodules  of  varying  size 
were  simulated  in  the  image  for  demonstration  purposes  before  the  noise  simulation  in  order  to  have  a  known  truth  for 
nodule  location  and  to  evaluate  the  effect  of  additive  noise  on  tumor  detection.  Noise  was  subsequently  added  in  order 
to  simulate  %,  14,  and  14  of  the  original  exposure  in  each  projection.  The  projections  were  then  processesed  with  the 
MITS  algorithm  to  produce  slice  images.  The  subjective  assessment  of  the  resulting  tomosynthesis  slice  images  show  a 
potential  decrease  in  dose  level  by  25-50%.  This  method  will  be  applied  to  a  study  of  dose  reduction  in  the  future  using 
human  subject  cases. 

Keywords:  digital  x-ray  imaging  (DX),  tomsynthesis  (RECON),  x-ray  tomosynthesis  (TSYN),  simulation  (SIM) 
Preferred  Presentation  Type:  Oral  Presentation 
Conference  Code:  MI03  -  Physics  of  Medical  Imaging. 


Introduction 

There  has  been  much  interest  in  dose  optimization  for  radiography  in  order  to  reduce  patient  exposure  while 
maintaining  adequate  image  quality  for  diagnosis.  However,  to  our  knowledge,  this  subject  has  not  yet  been  addressed 
for  digital  chest  tomosynthesis  imaging.  Much  of  the  misdiagnosis  in  projection  radiography  is  due  to  anatomic  noise 
from  overlying  structures.  Digital  tomosynthesis  is  a  method  which  reconstructs  longitudinal  planes  within  a  patient 
from  a  set  of  digital  projection  images  taken  over  a  limited  angle.  Tomosynthesis  increases  detectability  of 
abnormalities  by  removing  the  overlying  anatomy  and  improving  the  conspicuity  of  in-plane  structures.  Although 
tomosynthesis  has  been  around  for  many  decades,  it  has  only  recently  become  clinically  practical  with  the  advent  of 
flat-panel  detector  technology.  In  a  pilot  study  performed  in  our  lab,  50%  of  CT  confirmed  nodules  were  found  in  the 
PA  chest  radiograph  versus  81%  in  a  tomosynthesis  image  set2.  The  current  default  dose  setting  for  our  tomosynthesis 
imaging  system  is  roughly  equivalent  to  a  screen-film  lateral.  This  was  chosen  because  of  the  similarity  in  image 
quality  of  tomosynthesis  images  taken  at  this  exposure  to  conventional  chest  radiographs.  There  is  potential  to  lower 
tomosynthesis  exposures  because  the  conspicuity  of  nodular  opacities  in  the  tomosynthesis  images  at  the  current  default 
dose  setting  is  not  likely  limited  by  x-ray  quantum  noise.  Thus,  the  increased  detectability  of  abnormalities  provided  by 
tomosynthesis  may  decrease  the  need  for  the  default  dose  setting  which  gives  the  appearance  of  a  traditional  exposure. 
In  previous  studies,  the  addition  of  computer  simulated  noise  to  real  human  data  has  been  implemented  as  a  way  to 
investigate  the  diagnostic  accuracy  as  a  function  of  dose  reduction  ''".  The  methodology  discussed  in  this  work  provides 
a  realistic  simulation  of  lower  exposure  by  adding  stochastic  noise  that  has  been  filtered  by  the  characteristic  noise 
power  spectrum  (NPS)  of  the  system  to  the  projection  images  prior  to  tomosynthesis  reconstruction.  For  the  purpose  of 
this  study,  the  use  of  additional  simulated  lung  nodules  serves  to  provide  a  larger  sample  size  than  available  with  CT 
confirmed  nodules  as  well  as  the  ability  to  place  them  in  obscured  and  un-obscured  regions  to  in  order  to  evaluate  the 
effect  of  exposure  on  nodule  detection. 


Methods 

The  ultimate  purpose  of  this  methodology  is  to  simulate  an  image  acquired  at  a  reduced  dose  with  our  imaging  system. 
In  order  to  simulate  noise  on  previous  acquired  data,  which  realistically  emulates  the  actual  noise  of  the  system,  many 
steps  had  to  be  taken.  The  noise  power  spectrum  (NPS)  is  one  of  the  most  common  metrics  which  describe  the  noise 
properties  of  imaging  system.  We  experimentally  determined  the  NPS  of  the  tomosynthesis  acquisition  system  and 
utilize  it  to  filter  an  image  of  random  noise.  This  noise  is  added,  after  some  further  modifications,  to  the  original  image 
and  the  procedure  results  in  an  image  which  simulates  an  image  acquired  at  a  reduced  exposure. 


NPS  measurement: 

In  order  to  characterize  the  noise  properties  of  the  tomosynthesis  system,  the  noise  power  spectrum  (NPS)  was 
determined  experimentally.  In  the  spatial  frequency  domain,  the  NPS  is  the  variance  per  frequency  bin  of  a  stochastic 
signal.  In  order  to  correct  for  the  gain  of  the  system,  the  NPS  is  divided  by  the  square  of  the  mean  pixel  value  of  the  area 
under  analysis.  This  is  commonly  called  the  normalized  noise  power  spectrum  (NNPS)  )10'  n. 

NNPS(u,v)  =  NPS(u,v)  /  (large  area  signal)2  (1) 

A  previously  published  method  was  utilized  to  calculate  and  determine  the  NNPS10  n. 

At  each  dose  specified,  we  acquired  ten  tomosynthesis  image  sets  which  gave  us  10  flat  field  image  projections  for  each 
angle.  At  each  angle,  the  ten  flat  field  images  were  averaged  together  in  order  to  create  an  image  which  contained  the 
structural  noise  inherent  in  the  imaging  system.  This  averaged  image  set  was  subtracted  from  one  of  the  ten  original 
image  sets  in  order  to  remove  any  artifactual  background  shading  which  would  influence  the  NPS  measurement.  This 
background  subtracted  image  set  was  utilized  for  all  NPS  measurement.  In  order  to  account  for  the  alteration  in 
variance  due  to  the  image  subtraction,  the  resultant  NNPS  from  this  methodology  was  multiplied  by  N/N-l. 

A  640x640  pixel  area  in  the  center  of  each  background  subtracted  image,  segmented  into  128x128  overlapping  regions 
of  interest  (ROI)  was  used  for  analysis.  In  order  to  account  for  small  variations  in  regional  exposure  over  the  flat  field 
image,  each  ROI  was  scaled  by  a  ratio  of  its  mean  pixel  value  and  the  mean  value  of  a  designated  ROI. 

According  to  Saunders  et  al.10  the  NNPS  is  defined  according  to  the  following  relation: 


NNPS = - 
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Where  dA  is  the  pixel  area,  M  is  the  number  of  regions  used  for  analysis,  N  is  the  number  of  pixels  along  one  edge  of  an 
ROI,  ROI;  is  a  region  of  interest  in  the  area  under  analysis,  ROIi  is  the  ROI  in  the  top  left  corner  of  the  area  under 
analysis,  and  <ROI;>  was  the  mean  of  ROI;. 

Since  we  found  experimentally  that  the  NNPS  did  not  vary  greatly  between  the  different  projections,  we  averaged  all  of 
the  projections  in  order  to  give  the  final  NNPS  measurement  for  the  dose  under  analysis.  For  the  tomosynthesis  NNPS 
measurement,  the  NNPS  was  calculated  as  above  for  each  projection  image  and  the  average  NNPS  was  used  for  the  rest 
of  our  routine. 
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NNPS,  is  the  NNPS  from  a  single  projection  image,  and  K  is  the  number  of  projections. 


Noise  simulation 


In  order  to  create  additive  noise  which  has  the  same  noise  properties  as  the  true  projection  images,  the  texture  of  the 
noise  must  be  adjusted  according  to  the  experimentally  determined  NNPSt .  The  NNPSt,  from  only  one  dose  was  used 
to  filter  the  noise  for  all  dose  levels  because  the  shape  of  the  NNPST  was  found  to  be  relatively  consistent  at  the  range 
of  exposures  measured.  A  polynomial  curve  fit  was  applied  to  the  NNPSt  and  this  was  utilized  to  create  a  radially 
symmetric  estimated  NNPST  ( <NNPSt> ).  The  square  root  of  this  two  dimensional  <NNPST>  is  multiplied  by  the  two- 
dimensional  FFT  of  an  uncorrelated  Gaussian  noise  array  with  zero  mean  and  unit  variance.  The  noise  image  is  then 
converted  back  into  the  spatial  domain  by  taking  the  inverse  FFT.  The  resulting  filtered  noise  now  has  frequency 
content  consistent  with  the  tomosynthesis  imaging  system.  Before  adding  the  filtered  noise  image  to  the  original  image 
(I original ),  the  variance  of  the  noise  has  to  be  adjusted  so  that  the  resultant  image  ( I  simulated )  has  a  signal  to  noise  ratio 
(SNR)  consistent  with  a  reduced  dose  image. 

As  radiation  dose  decreases,  the  SNR  of  the  image  also  decreases.  However,  since  we  are  utilizing  previously  acquired 
clinical  data,  there  is  no  way  to  decrease  the  dose  experimentally  in  order  to  find  the  SNR  of  a  reduced  dose  image. 
Therefore,  we  added  an  image  containing  filtered  noise  (/„„,„,)  to  the  original  projection  image  ( hnginai )  that  decreased 
the  SNR  appropriately  to  achieve  the  target  reduced  SNR  (SNRreduced).  The  result  is  an  image  which  simulates  a  reduced 
exposure  (Isimuiated)- 
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simulated 


I  original  ^  noise 


(3) 


The  target  SNR  cannot  be  achieved  without  considering  the  noise  already  present  in  the  original  image  (/original)- 
Therefore,  because  variances  add  for  sums  of  uncorrelated  random  variables,  the  required  variance  to  simulate  a  lower 
exposure  can  be  calculated. 


^ simulated  ^ original  ^ noise 

The  relationship  between  variance  in  the  image  (o2)  and  the  mean  pixel  intensity  (<  I  >)  was  experimentally  determined 
in  order  to  account  for  the  spatial  variation  in  photon  flux.  From  this  relationship,  we  derived  that  SNR  varies  as  a 
function  of  exposure  (snr(E)).  This  was  done  assuming  that  the  expected  pixel  value  of  the  image  (<  /  >)  is 
proportional  to  exposure  in  a  linear,  quantum-limited  detector. 

SNRoriginal  =  =  snr(Eoriginal ) 
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SNRreduced  =  Snr(E  reduced  ) 
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Because  we  are  adding  noise  to  the  original  image  in  order  to  attain  a  reduced  SNR,  we  can  solve  for  the  additive  noise 
variance. 
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The  filtered  noise  is  multiplied  by  the  variance  determined  from  Equation  8.  This  process  adjusts  Inoise  for  a  spatial 
variance  consistent  with  the  local  pixel  intensity  in  the  original  image  in  order  to  account  for  varying  attenuation  from 
different  anatomical  structures.  This  process  results  in  a  final  noise  image  Inoise  which  contains  the  same  frequency 
content  as  the  inherent  system  noise  and  is  also  adjusted  for  the  local  spatial  variations  in  exposure.  This  final  Inoise  is 
added  to  Ioriginai  to  create  the  reduced  exposure  projection  image  Isimuiated  which  has  and  SNR  equal  to  SNRreduced. 

Nodule  Simulation: 

In  each  set  of  projections,  five  subtle  lung  nodules  ranging  from  4  to  8  mm  in  diameter  were  simulated  at  different 
locations  in  the  left  lung  on  the  original  projection  images12' 13.  The  locations  of  the  nodules  were  chosen  carefully  to 
evaluate  the  effect  of  exposure  level  on  the  detectability  of  nodules  under  various  anatomical  conditions:  mediastinum, 
vessel  crossings,  and  behind  ribs.  The  location  of  each  simulated  nodule  on  each  projection  imagewas  calculated  using 
the  fractional  angular  displacement  from  horizontal  (0°  is  located  at  projection  36).  The  resulting  projections  with  five 
simulated  nodules  were  then  passed  into  the  aforementioned  noise  simulation  routine.  Nodule  simulation  is  for 
demonstration  purposes  only.  It  is  utilized  to  have  a  known  truth  of  tumor  location  and  to  evaluate  the  effect  of  dose 
reduction  on  nodule  detection;  but  is  not  necessary  for  the  typical  application  of  this  methodology  to  assess  the  clinical 
detection  performance  of  reduced  dose  tomosynthesis  images. 

Acquisition  method 

A  prototype  system  constructed  in  our  laboratory  with  a  commercial-grade  41x41  cm  a:Si/Csl  flat-panel  radiographic 
detector(GE  Healthcare,  Milwaukee,  WI)  and  rapid  acquisition  hardware  was  utilized  for  all  measurements  and  image 
acquisition.  NPS  measurement  was  done  in  accordance  with  the  International  Electrotechnical  Commission  (IEC) 
standard  of  measurement  (IEC  61267  Ed.  2.0  b:2005;  RQA9  technique)9  slightly  adjusted  for  tomosynthesis  image 
acquisition.  The  antiscatter  grid,  system  faceplate,  and  0.2  mm  of  Cu  filtration  were  kept  in  place  for  system  calibration 
and  NPS  measurement.  Beam  quality  was  achieved  with  120  kV  and  39.9  mm  A1  filtration  (half-value  layer  of  11.5 
mm).  All  flat  field  image  sets  and  exposure  measurements  were  made  with  the  standard  tomosynthesis  acquisition 
procedure  which  utilizes  71  projection  images  over  20  degrees  and  acquired  at  6.4  frames  per  second  with  the 
manufacturer  supplied  gain  and  offset  corrections  applied.  All  images  were  acquired  with  39.9  mm  of  A1  and  0.2mm  of 
Cu  filtration  in  place.  Measurements  were  made  at  four  exposure  levels  using:  0.32,  0.4,  0.64,  and  1.25  mAs/projection. 

Due  to  the  geometry  of  the  system,  it  was  not  possible  to  place  the  ionization  chamber  (MDH  Model  1015,  10X5-6 
ionization  chamber,  Radcal,  Monrovia,  CA)  in  the  beam  during  acquisition  of  the  NPS  images.  Therefore,  the 
cumulative  exposure  at  each  dose  level  for  all  71  projections  was  measured  at  the  plane  of  the  detector  utilizing  a  small 
field  of  view,  after  the  NPS  image  acquisitions  for  each  dose  were  completed. 

In  order  to  give  examples  of  how  this  dose  reduction  simulation  method  performs,  tomosynthesis  datasets  were  acquired 
of  real  human  subjects  at  120  kVp  and  an  mAs  setting  which  was  determined  based  upon  either  subject  thickness  or  a 
photo-timed  value  using  the  Automatic-Exposure  Control  (AEC)  from  a  clinical  digital  chest  radiograph  unit  (GE 
Healthcare,  Milwaukee  WI).  The  range  of  mAs  values  was  0.32  to  1  mAs  per  projection. 

After  nodule  and  noise  simulation,  69  planes  were  reconstructed  from  each  set  of  the  projections  with  simulated  noise 
and  nodules  using  the  MITS  (Matrix  Inversion  Tomosynthesis)  algorithm  l4' 16  developed  in  our  lab,  with  5  mm  plane 
spacing.  A  sliding  average  of  7  adjacent  planes  was  used  to  improve  image  noise  and  resulted  in  the  final  image  set 
which  was  utilized  for  the  analysis  of  the  simulated  dose  reduction  method. 


Results 

The  normalized  NPS  (NNPS)  from  four  different  dose  levels,  which  are  representative  of  what  are  used  during  clinical 
acquisition,  were  analyzed:  0.32,  0.4,  0.64,  and  1.25  mAs  per  projection,  with  the  cumulative  exposure  from  a  full 
tomosynthesis  acquisition  of:  1.83,  2.35,  3.85,  and  7.8  mR  respectively.  Figure  la.  shows  that  the  shape  NNPS  did  not 
vary  greatly  at  different  projections  and  therefore  we  averaged  them  to  form  the  composite  NNPS  ( NNPSt ,)  for  each 
dose.  A  fifth  order  polynomial  curve  fit  was  utilized  to  estimate  an  equation  for  the  NNPSt.  The  resultant  fit  is  shown  in 
Figure  lb.  along  with  the  NNPST.  Figure  lc  shows  the  NNPS  from  one  projection  angle  (0°)  for  each  of  the  dose  levels. 
Figure  Id  demonstrates  that  the  system  has  a  very  low  level  of  electronic  noise  because  the  product  of  NNPS  and 
exposure  does  not  vary  greatly  over  the  range  of  exposures  evaluated. 


NNPS  For  one  Dose  Level  NNPS  For  One  Dose  Level 


Figure  1:  (a)  Ensemble  of  normalized  NPS  (NNPS)  curves  from  one  dose  acquisition  of  all  71 
projection  images,  (b)  The  averaged  NNPS  curve  from  one  dose  acquisition  with  the  associated 
curve  fit.  (c)  Example  NNPS  of  all  four  exposures,  (d)  The  NNPS  multiplied  by  the  exposure  level 
(curves  would  all  be  the  same  if  the  detector  was  quantum-limited). 

The  linearity  of  the  detector  was  found  by  computing  the  spatial  mean  and  variance  of  the  flat  field  projections  utilizing 
the  central  80%  of  the  image  for  each  exposure  level.  (Note:  Variance  was  determined  from  flat  field  projections  with 
the  background  subtracted  and  corrected  with  N/N-l).  Using  this  relationship,  the  amount  of  noise  in  the  reduced  dose 
image  is  computed  for  each  pixel.  In  our  methodology,  we  utilize  the  characteristic  response  curves  of  the  system 
shown  in  Figure  2a  in  order  to  estimate  the  relative  exposure  level  in  each  pixel  of  the  clinical  image  set.  From  this 
relationship,  we  utilize  the  relationship  between  the  SNR  and  exposure  level  and  can  determine  the  SNR  at  the  target 
dose  level.  Figure  2  (a,c,e)  demonstrate  these  empirically  derived  relationships.  From  Figure  2  (b,d)  it  can  be  seen  that 
these  relationships  are  not  highly  dependent  on  projection  angle.  We  averaged  the  mean  and  variance  over  all 
projection  angles  in  order  to  derive  the  relationships  for  the  system  performance  due  to  exposure. 
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Figure  2:  (a)  Response  of  acquisition  system  relating  mean  pixel  value  (MPV)  and  exposure,  (b)  MPV 
as  a  function  of  projection  angle,  (c)  Empirical  relationship  between  variance  and  exposure. 

(d).  Variance  as  a  function  of  projection  angle,  (e)  SNR  as  a  function  of  exposure. 


Figure  3:  Simulated  nodules  which  are  added  to  the  projection 
image  prior  to  noise  simulation  and  reconstruction 


(a) 


(c) 

Figure  4:  (a)  Full  dose  projection  image  at  0°.  (b)  Simulated  half  dose  projection  image  at  0°.  (c)  Full  dose 
reconstructed  tomosynthesis  image  (planel4).  (d)  Simulated  half  dose  reconstructed  tomosynthesis  image 
(plane!4). 


Nodule  simulation  was  performed  on  each  projection  image  prior  to  noise  simulation.  An  example  of  the  simulated 
nodules  that  were  added  to  the  original  image  is  shown  in  Figure  3. 


The  resultant  projection  images,  after  the  noise  simulation  method  was  applied,  approximate  the  correct  SNR  for  an 
image  taken  at  a  reduced  exposure.  This  noise  propagates  through  the  tomosynthesis  reconstruction  algorithm  and 
simulates  how  reduced  dose  would  affect  nodule  conspicuity  and  detection  accuracy.  Figure  4  shows  the  effect  of 
adding  noise  to  simulate  half  the  exposure  of  the  original  image  in  both  the  original  projections  and  the  resultant 
tomosynthesis  reconstructions. 

Figure  5  demonstrates  the  nodule  conspicuity  dose  response  under  different  conditions.  Each  ROI  shown  is  16  cm2  with 
the  simulated  nodule  located  in  the  center.  The  same  ROI  location  was  used  for  each  row.  The  first  row  shows  a  4  mm 
simulated  nodule  placed  in  the  lower  left  lung;  the  second  row  shows  a  6  mm  simulated  nodule  placed  behind  a  rib;  and 
the  third  row  shows  an  8  mm  simulated  nodule  placed  in  the  mediastinum. 


Figure  5:  The  effect  of  low  dose  simulations  on  nodule  consipicuity. 

1st  column  shows  the  ROI  from  the  projection  image  at  100%  dose; 

2nd  column  is  the  ROI  from  the  tomosynthesis  reconstructed  plane  at  100%  dose; 
3rd  column  is  the  ROI  from  the  tomosynthesis  reconstructed  plane  at  75%  dose; 
4th  column  is  the  ROI  from  the  tomosynthesis  reconstructed  plane  at  50%  dose; 
5th  column  is  the  ROI  from  the  tomosynthesis  reconstructed  plane  at  25%  dose. 


Initial  subjective  analysis  of  the  simulated  reduced  dose  images  suggest  that  this  method  adequately  emulates  a 
tomosynthesis  image  acquired  at  a  reduced  dose.  This  method  will  be  used  in  a  future  study  to  determine  the  optimum 
reduced  exposure  for  tomosynthesis  image  acquisition. 


Discussion 

In  this  study,  we  made  the  assumption  that  the  shape  of  the  NNPS  would  be  independent  of  exposure.  There  is, 
nevertheless,  some  variation  due  to  exposure  and  an  additional  step  can  be  implemented  which  may  slightly  improve  the 
performance  of  the  method  presented  in  this  work.  However,  with  the  small  amount  of  variation  at  these  low  exposure 
levels,  only  a  minor  improvement  can  be  expected. 

Conclusions: 

The  methodology  described,  generates  tomosynthesis  images  subjectively  equivalent  to  images  acquired  at  reduced 
dose.  When  coupled  with  simulated  nodules,  this  approach  may  be  used  with  human  observers  to  conduct  ROC  studies 
of  observer  performance  for  nodules  in  both  obscured  and  un-obscured  lung. 

Future  Work: 

Clinical  subject  data  with  simulated  noise  and  nodules  will  be  used  in  an  ROC  study  with  chest  radiologists  in  order  to 
more  accurately  determine  the  detection  accuracy  in  tomosynthesis  images  at  reduced  dose  levels  for  dose  optimization. 
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